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O ' 

O ' Abstract 

, We consider the singular configurations of gravitating gas [IJ wliich 

Q I could be used as a model for disk galaxies. The simplest steady con- 

■ figuration, which corresponds to rotation of stars around center gives 

flat rotational curve, provided the density distrubution has a tail, de- 
caying as 1/r. This result based exclusively on Newtonian gravity and 
does not involve dark matter. 



> 

m 

m ■ 1 Introduction 



In the present paper we continue started in [T] study of the dynamics of the 
O ■ system called gravitating gas. This system is a kind of continuous media 

which is constituted by the particles (stars) interacting with each other only 
gravitationaly i.e. according Newton's law. As a mechanical system it could 



Oh' 



O I be considered as a continuous limit of the A^-particles system described by 

' Hamiltonian: 

^ H = j2^'^ J2 ^^'"^^ (1) 

_c3 \ where are canonical coordinates of particles (stars) with masses mj, 

G-gravitational constant. 

Speaking about galaxy, where ~ 10^^ it seems reasonable to consider 
for the system ([T]) the limit N ^ oo (with assumption = m) and substi- 
tute the configuration of particles by continuous distribution of its density 
and velocity, as it is done in the theory of gas or fluid. 

Certainly, we have to have in mind that the resulting system- gravitating 
gas is a very strange gas due to attractive interaction of its constituents. In 
the usual gas or fluid, the constituents (atoms or moleculars) have a repul- 
sive interaction and being put in a volume it spread in space, filling after 
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some time the whole volume with uniform density. As it was proved in 
[1], the gravitating gas, because of attractive interaction may form isolated 
steady configuration with asymptotically vanishing density. In particularly, 
the galaxy named Hoag's object may be an example of such gravitational 
soliton 

Here we are going to consider another kind of steady configurations of 
gravitating gas which has a singular distribution of density 

p{x) = 5{z)p{r,<P), (2) 

where z, r, are cylindrical coordinates. This kind of solutions may be con- 
sidered as an idealized model of disk galaxies. 

Apart from the application on the galactic scale these singular solutions 
may be used for the explanation of planetary rings such as Saturn ring with 
its surprisingly thin structure formed exclusively by gravity. 

2 Equations of motion 

In fluid (gas) mechanics there are two different pictures of description. The 
first, usually refereed as Eulerian, uses as the coordinates the space dependent 
fields of velocity and density. The second, Lagrangian description, uses the 
coordinates of the particles x{^i,t) labeled by the set of the parameters 
(the numbers of particles in ([T])), which could be considered as the initial 
positions ^ = x{^i,t = 0) and time t. The useful physical assumption is 
that the functions x{C,i, t) define a diffeomorphism oi D <0 and the inverse 
functions ^{xi,t) should also exist. 



•^j (^« 5 ^) 



^j{xi,t) ^ ^, = ^j- (3) 
The density of the particles in space at time t is 

p{x,t)= [ £^po{Ci)5{x- x{^i,t)), (4) 



where po(0 is the initial density at time t = 0. The velocity field v as a. 
function of coordinates x and t is: 

v{xi,t) = x{^{xi,t),t), (5) 
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where ^{x,t) is the inverse function ([3]). The velocity also could be written 
in the following form: 

^/ ,^ ^ / d^^Po{^i)xi^u t)Six - Xj^i, t)) 

or 



p{xi, t)v{xi, t) = J d ^po{^t)x{^i, t)6{x - x{^i, t)). (7) 
Let us calculate the time derivative of the density using its definition (jl]) : 

/O 

= -i / d'^Po{^^m^i, t)S{x - f t)) 

d 

= -—p{xi,t)v{xi,t) (8) 

In such a way we verify the continuity equation of fluid dynamics: 

p{xi, t) + d{j){xi, t)v{xi, t)) = 0. (9) 

Using the coordinates x{C,i,t) as a configurational space variables we can 
write the Lagrangian for the continuous generalization of the system, de- 
scribed by ([ID 

,3. ^ mx\^i,t) Gm^ r pQ^ji) Po{i'i) .^ 

d epote) — ^ — + ^jdid^ \m^.t)-mit)y 



n.m,t) + On.^jd^i'p.^O^§§^§;^, ^ 0. (11) 



The equations of motion, which follow from the Lagrangian (ITOl) have the 
form: 

x{^i,t) -x{^'i,t) 

Now we need to translate the equations ( ITTi) on to the language of Euler 
variables. For this let us differentiate both sides of the equation(l6]) with 
respect to time: 

d f " 

—p{xi,t)v{xi,t) = J d^^pQ{^i)x(^i,t)S{x - x{^i,t)) 

+ 1 d'^poi^.mu t)^^S{x - f te, t)). (12) 



Substituting x{^i,t) from the equation (ITT]) and transforming the second 
term, as we did in derivation of the equation ([7])we obtain: 

d d 



dt 



dxk 



d'iPo{ii) -Gm / d^e'Po(e 



6{x-x{^i,t)) 



(13) 



To transform the last integral we first perform integration over the ^: 



d'M^i) -Gm / d'^'poii 



-Gmp{x„t) I d'CpoiCi 



6{x-x{Ci,t)) 



\x-m,tw 

Now let us insert in the integral over ^' the unity 

d'ySiy-x{C,t)) 
and change the order of integration. Finally we arrive at 



(14) 



(15) 



d'M^^) -Gm / d'ePo{e^ 



^1 /WN X{ii,t) - X{^'i,t) 



\x{^,,t)-x{^'„t)\^ 



5(f-x(C',t) 



d 



Gmp{xi,t)— I d y 



p{yi,t) 

\x-yl 



(16) 



The equation of motion ffTTl) in terms of Euler variables takes the following 
form: 



d d 

p{Xi, t)v{Xi, t) + (p{Xi, t)v{Xi, t)Vk{Xi, t) 



dt 



dxk 

Gmp{x,,t)^ f d'y^^. 

ox J \x — y\ 



(17) 



Using the continuity equation we can rewrite this equation in more simple 
form: 



But this is possible only if density does not vanish in some domain. If for 
example we want to study configurations with density like given by ([2]), we 
must keep the equation of motion in the form f|T7|) . 

Indeed, let us consider the configuration of our system with velocity and 
density of the type 

p{x,t) = 5{z)p{x,y,t) 
Vi = {vi{x, y, t),V2ix, y, t), 0). (19) 

Substituting (fT9l) into (fTTI) we obtain: 
d 

—S{z)p{x, y, t){vi{x, y, t),V2{x, y, t), 0) + 

d 

^6{z)p{x, y, t)vaix, y, t){vi{x, y, t),V2ix, y, t), 0) 



dx^ 

G m S(r.)p{x, y, t)(^, 1. 1) / d.^^— ^^^^L— . (20) 

ux uy uz J ^ ) + [y ~~ y j + -2 

where index a takes the values 1,2. Apparently, the derivative of integral 
with in the r.h.s of (l20l) with respect to z vanishes after multiplication on 
S{z), because the integral is even function of z, so the third component in 
both sides of (!20|) disappears and after making use of the continuity equation 
we obtain 

— t;a(x, y, t) + vi,[x, y, ^) = 

A / dx'dy' ^ ^("'^'^^ (21) 



dx„ 



x — x'y + {y — y'Y 



where indexes a, h take the values 1,2. This equation has to be completed 
with continuity equation, which now has the following form: 

d d 

— p(x, y^t) + — (p(x, y, t)va{x, y, t)) = (22) 

Note, that as it should be anticipated, the equations fl2Tl) and fl22|) took the 
usual form of 2-dimensional gas dynamics, but the potential inherited its 
3-dimensional form. 
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3 Simplest steady solution and its rotational 
curves. 

Here we are going to consider the time independent configurations of gravi- 
tating gas, tlierefore all time derivatives in the equations fl^Tl) and fl2^ will 
disappear. 

Vb{x,y)—Va{x,y) = Gm- — / dx'dy'—^ ' , 
dxh dxaJ ^{^x-x'Y + {y-y'Y 

d 

— (p(x, y, t)va{x, y, t)) = 0. (23) 

One of the simplest configuration of gravitating gas describes the rotation of 
the particles (stars) around central point. For this configuration we have: 



p = p{r), r = \Jx'^ + y'^. (24) 

This ansatz satisfies the continuity equation and the main equation give a 
relation between velocity v{r) and density p(r): 

v'^{r-) = -GMrdr J dr'r' p{r')t{r,r'), (25) 

where the kernel t{r,r') is given by: 

t{r,r')= dcj) , K(A;), (26) 

Jo yr^ + r'^ — 2rr'cos(j) r + r' 

k 



(r + r'Y 

and K(/c) the complete elliptic integral of the 1-st kind. So, from this solution 
follows that any distribution of density is possible and and we have to consider 
the stability of solution in order choose p{r). 

The equation (l25l) shows that is in case of flat disk configuration of stars 
the famous Newton's theorem is not valid. Of course, this theorem which 
states that |[2j 

v^{r) = Gm^^ (27) 



where M(r) is the total mass inside the orbit, is true only for circular orbits 
in spherically symmetric density distribution. In our case we can rewrite (1251) 
in the following form using modular transformation of modulus k: 

. rr A r' 4 r \ 

v^(r) = -Gmrdri dr'r' p(r')-K(-) + dr'r' p(r')-K(-)) . (28) 
\Jq r r Jr ^1 ^1 J 

The first term in the r.h.s. has the factor i, but dependence of r is present 
also in the integrands of both terms (the dependence from the integration 
limits is cancelled of course). 

So, what we can say about typical rotational curves i.e. the functions v{r) 
for different densities? Apparently, if the density is strictly localized within 
the circle with radius i?, when, for r 3> -R we obtain the usual result: 

t;2(r) = Gmy , (29) 

where M is the total mass (which all is in the circle r = R). 

In the disk galaxies almost all mass is concentrated near the center de- 
caying exponentially e~^/'^ with /i ~ 1 — 5 kpc [2]. In this case (l29ll could 
be applied for r h. But, in the same time the disk of the galaxy, as it seen 
is much bigger then h, therefore we may suggest that there is a tail in the 
density distribution which drops slowly, e.g. as a power of r. 

In order to explore the influence of such kind of densities, we need to 
know properties of the kernel t{r,r') in the integral fl25|) . Apparently it is a 
symmetric function of r, r' and has the logarithmic singularity when r —>■ r' 
(the singularity of elliptic integral K(fc) for k 1). Unfortunately there is 
only one case when we can express the integral fl2S]) in terms of elementary 
functions: for density 

p(r) = - (30) 
r 

the integral is given by 

J dr'r'p{r')t{r, r') = -c 27r In r + C, (31) 

where C is a divergent constant, which is inessential for calculation of v{r): 

v^{r) = G mc2n (32) 

So, the rotational curve for p(r), produced by the density fl30|) is flat! Cer- 
tainly, in order to have fiat rotational curve, it is sufficient to have the be- 
haviour of density given by (l30i) only asymptotically, starting e.g. from the 
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distances ~ 6 — 8kpc. For smaller distance the exponential distribution pre- 
vails, not changing the behavior at large distance. 

As it follows from the analysis of the integral in fl25l) . the density fl30l) 
is the limiting case for the convergence of integral which defines f^(r). If 
asymptotically 

P(0 = " > (33) 

then we have decreasing rotational curve. 



4 Discussion. 

The result we obtained in the previous section is not very unexpected and is 
the direct consequence of the flat density distribution. In order to explain it 
let us derive it in a more traditional way without referring to gravitational 
gas, considering the motion of particle in axially symmetric density distri- 
bution p{xi) = p{r, cos 9), where r,9,(l) now are spherical coordinates. Then 
the equation of motion is given by 

p(r', cos^^') 
i/r^ + r'^ — 2rr' cos B ' 

where is the angle between x and x': 



S{t) = Gmd f dr'r'^d(P'd cos 9' p{r', cos 9') ^ ^^^^ 
J \/r^ + r'^ — 2rr' cos 



cos = COS 9 cos 9' + sin 9 sin 9' cos(0 — (f)') (35) 

Now let us assume that the test particle moves along circular orbit with 
angular velocity uj{r), therefore 

f(t) = -uj\r)x{t) = -!^f(t), (36) 
from where we obtain: 

v\r) = Gmrdr I dr'r'^d<f)'d cos 9' , Pi^'^^^^^') ^ (37) 

J V + — 2rr' cos 

Using well-known decomposition of the square root in ( 134|) we can expand 
the r.h.s. of (1371) in the following form: 

f^(r) = Gmrdr ^ / dr' r''^ dcp' d cos 9' p{r' , cos 9')—^ Prt{cosQ) 

roc j^'^ -, 

+ / dr' r'^d(f)'d COS 9' p(r', COS 9')— —PnicosQ) . (38) 

Jr j-ln+l J 
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If the density is spherically symmetric, the integration over angles 9' and 0' 
leaves only the first term in the sum and we immediately recover Newton's 
theorem. In any general case, all terms of the sum are non-zero and we 
have complicated behaviour of rotational curve. If for example the density 
has the form we had considered above, i.e. p ~ p{r)8{rcos{6)), and the 
orbit lies in the plane xy we are coming to the equation from the previous 
section. Certainly, this derivation can not substitute the whole dynamics of 
gravitating gas we had considered, because the potential we used there is 
created by the same particles, which move in the potential. 

Now we can estimate the order of density we need to fit the observation 
data of the velocities on the fiat part of the rotational curves. For typical disk 
galaxy these velocities reach the 150 — 300A;m./sec. If we shall take as mass 
of the particle to be equal to the mass of our Sun — 2 ■ 10^" kg, the coefficient 
GM2t[ will be ~ 8,4 ■ 10^°m^/sec^. From here we obtain for the constant c 
the values (0, 27 — 1, 1) ■ 10^^°m~^, which corresponds to (0, 8 — 3, 3) ■ lO^pc"^ 
The surface density therefore becomes 



where r is measured in pc. On the distance of e.g. 10 kpc, this formula 
gives us the surface density ~ (0, 8 — 3, 3) • 10^ stars on square parsec. These 
numbers is not in contradiction with observational results 

The precise measurement of velocities of stars in various galaxies which 
was done during last 30 years [3],|1] made it possible to discover the peculiar 
behaviour of the rotational curves. This in turn led one part of community of 
astronomers and astrophysicists to the hypothesis of dark matter, the other 
to modification of gravitational law Later another confirmation of the 

existance of dark matter has emerged: the data on gravitational lensing and 
cosmic microwave background radiation. However, the argument, given by 
the rotational curves stays aside as a completely model independent measure- 
ment, not based on any assumptions of global structure of Universe. That is 
the reason why it is so important to explore all possibility within traditional 
physics and I hope that the present paper is one more step in this direction. 
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